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Abstract. Interacting systems consisting of two rotators and 
a point mass near a hyperbolic fixed point are considered, in 
a case in which the uncoupled systems have three very differ- 
■ ent characteristic time scales. The abundance of quasi periodic 

t-H [ motions in phase space is studied via the Hamilton-Jacobi equa- 

tion. The main result, a high density theorem of invariant tori, 
is derived by the classical canonical transformation method ex- 
tending previous results. As an application the existence of long 
\ heteroclinic chains (and of Arnol'd diffusion) is proved for sys- 

tems interacting through a trigonometric polynomial in the angle 
variables. 
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o : 

^ . §1. The system. 

oo ! 

I Let (01,0:2) — a G T 2 be a pair of angles; let A — (Ai,A 2 ) G R 2 be their 

"^5 ■ conjugate momenta or "actions", and let (p, q) G R 2 be a further pair of canonically 

conjugate coordinates. We consider the Hamiltonian function, depending on two 
*^ ■ dimcnsionless parameters, e, 77 > 0, defined by: 

U = fcfa'Ai) +t)-^oj 2 A 2 + G(pq,n?Ax) + e/fa'Ai, a,p,q) (1.1) 

We suppose that h, G, f are real analytic functions of their arguments in the domain 
fj'l^il < R, \p\, \q\ < r, a G T 2 for some R, r > 0. We also suppose that / is a 
, trigonometric polynomial of degree N in the variables a . 

For s = the Hamiltonian (1.1) represents, physically, a system consisting of two 
rotators and a point mass near an unstable fixed point (p = q = 0). 



The instability time scale of the unstable equilibrium is g = g(pq,rji A\) = 
d x G(x, n? Ai)\ x=pq — 0(1) if, as we shall suppose, g > 0. 

Define rjiuji = ?p 0Ji(pq, ri^Ax) = c^i [M^^i) + G(pq,rji Ai)]: the two rotators 
rotate with angular velocities ij'wi and n~^uj2 with ratio 0{n) if, as we shall suppose, 
< Lo a < |wi I < &b for suitable constants U) a ,u)b. 

Hence the system has three time scales of respective orders n~i , 1, called respec- 
tively: slow (rotator #1), normal (unstable motion), fast (rotator #2). We suppose 
i] < 1 so that the fast rotator really is the isochronous one and the slow rotator is 
the anisochronous one. 

§2. Characteristic parameters. 

A few characteristic parameters can be associated with the system because of the 
analyticity assumptions. 

Let 2po, 2£o, (2kq) 5 ) (2kq) 5 , 2kq measure the size in the complex planes of a holo- 
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morphy domain in the Ai,cti,p,q and x — pq variables, respectively. We denote: 

Q p ={|ReA!| < v -i(R + p), |ImAi|<»j-ip} 

Ut={\lm aj \ <£}, ^ = {N,|g| S K = {M< K } 

Our analyticity hypotheses imply the existence of po,£o, K o > such that h(r)iAi) 
is holomorphic in Q,2p > 772^1) is holomorphic in 52« x Q 2 p an( l / m Q2p x 

t^2Co x ViKo = ^2p ,2{o,2K - 

Suppose: 

(i) .g( P (Z,r/i J 4 1 ) d i: / a x G(x,ry*A 1 )| ;c=P9 ^ 0; 

(ii) u\(pq, rjiAi) d = dh(r)?Ai) + dG(pq, r]^Ai) ^ 0, where 9 denotes differentiation 
with respect to the argument a = 772 A\\ 

(iii) mi{pq, Tj^Ai) d = d 2 h(r]iAi) + d 2 G{pq, rj^Ai) ^ 0. 
Hence we can define To, Eq, Mq > so that: 



< A =min{p ,Ko}, |/|, \h\, \G\ < E 

< r < M M , 77^)1, |w 2 |, |fl(pg, 77*Ai)| < SqAq 1 (2.2) 
< Mo" 1 < |mi(pg,»?iAi)| 

in the holomorphy domain T>2 POl 2£, ,2n ■ We suppose rj < 1 (see §1) and, for simplicity, 
£o,£ < 1. 

The parameters So, To, Mo, Ao, £0 will be called characteristic parameters of the sys- 
tem (note that roAo,WbAo < Eo). Only the dimensionlcss parameters £0, EoTq Xq 
and Mo-EoAq 2 will control the constants in the main bounds (see (4.3) and (5.5)). 

§3. Hamilton Jacobi equation. 

The Hamilton- Jacobi equation for the Hamiltonian (1.1) will be an equation with 
two unknown functions $(772^, a,p',q) and G(x, 772 A[) given by: 

+ ai $)) + v'^MA + d a2 f>) + G((p' + d q $)q, r}i(A[ + ai $)) 
+ e/^* (Ai + ai $), a,p' + d q <S>, q) = (3.1) 
= h(^A[) + rj" W 2 + G(p'(q + d p ,$),^A[) + G(p'(q + c^),tM) 



If the functions $, G existed one could set: 

n = n' 4- f>,.& 

(3.2) 



A=A' + da$, p = P ' + d q <P 



a'=a+dA'®, q' = q + d p '& 

and the change of variables (A, a,p,q) < — > (A', a',p',q'), where defined, would 
cast the Hamiltonian (1.1) into the form: 

Ho = + rj" W 2 + G(p'q' ,r,i A[) + G(p'q', rM) (3.3) 

thereby integrating the system, because the motions in the new coordinates would 
become A' = const., a ' — > a' + w 't, p' — > p'e~^ * and q' — ► q'e^ * with: 

5 / = 0V,^^i)=^G(a ; / ,7 ? 5Ai) + ^G( a; ',7 ? 5A' 1 ) (3.4) 
w' = u'(x?,ti*A[) = (^^^(ry^Ai) +7 7 ^'G(x',773A 1 ) +772a'G(x',775A' 1 ),7 ? -^ 2 ) 
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where x' = p'q' and & denotes differentiation with respect to the argument a' = 

This is in fact, as it is well known, not possible in general. However one proves the 
following result. 

Proposition. There exists a canonical change of variables like (3.2) such that 

(a) it is of class C°° in all variables and in e for e small enough; 

(b) it is close to the identity within 0{e) in C°°; 

(c) there is a function G(x' , A[ ) which is of class C°° in all variables, including e, 
which is divisible by e and such that in the new variables the Hamiltonian coincides 
together with its derivatives, with (3.3), on the set (without interior points) of the 
(Ai,p,'q', a')'s such that the rotation vector lu' = tj'(p'q' ,rp A[) defined in (3.4) 
verifies the Diophantine condition: 

k' • E\ > C{n)\v\- T for all 0^ v_e7? (3.5) 

with t > 1 prefixed and C{n) > 0, provided e is small enough, depending on C(r)) 
and the other parameters of the model. 

The above proposition is essentially proven in [CG], §5, Lemma 1', (see [G2] and [P] 
for cases without gaps). The proof in [CG] deals really with the subset of the set the 
variables A 1 and x' — p'q' where u '(p'q 1 , rp A[) = (1 + ^(p'q 1 , r]^ A[)) a>'(0,rp A[) 
for some "f(x', a') of class C°° in its arguments: by examining the proof one sees 
that the result holds under the above more general condition. 

Remarks. (1) Calling W — W et c(rf\ the set in the phase space such that the corre- 
sponding rotation vectors u±' verify (3.5), one says that the Hamilton- J acobi equa- 
tion is soluble on W and, in W, casts the Hamiltonian (1.1) into the normal form 
given by the r.h.s. of (3.3). 

(2) Since the function oj' has non zero gradient with respect to A' (under the 
mentioned condition of smallncss of e; sec (2.2)), the volume W of phase space 
where (3.5) holds has a complement with measure bounded proportionally to C{rf). 
And at fixed x' = p'q' the measure of the subset of the interval [— n~ 2 R i -q- 2 R\ on 
the A[ axis has also measure bounded proportionally to C(n). The latter subset 
consists of a sequence of small intervals (whose total length has size < 0(C(n))) 
that are called "gaps" for natural reasons. 

(3) The values of A[ such that (3.5) holds with x' = p'q' = are points where the 
Hamiltonian coincides with one like (3.3) together with its first derivatives. Thus 
they clearly correspond to invariant tori (if p' = q' = 0, also called "persistent 
tori") and to their stable and unstable manifolds (if q' = or p' = 0, respectively, 
sometimes called whiskers). On such manifolds the motion takes place following the 
evolution described after (3.3). 

If we simply apply the results in [CG], the bounds on the gaps between the invariant 
tori (which persist notwithstanding the perturbation) are found to be "too large" 
with respect to the li homoclinic splitting" size (it has become usual to say that the 
splitting is "exponentially small in n~ 2" as 77 — > 0). The latter, see below (remark 
3 in §5), is the quantity of main interest in the theory of drift and diffusion. A 
simple way to use this result to discuss diffusion in phase space is to choose e very 
small: but it turns out that it has to be chosen too small for the purpose of several 
interesting applications (and one says that it has to be chosen "exponentially small 
with n~ 2"); we postpone to §5 a more technical discussion. 

This paper is dedicated to overcome the limitations that have been briefly alluded 
to and to show that the restrictions to be put on the size of e can be greatly improved 
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with respect to the ones that follow directly from previous work, (§5 and Appendix 
A10 of [CG]). Therefore we shall proceed in a somewhat different way with respect 
to the one outlined above, and we shall perform a preliminary step (§4) before 
applying the analysis of [CG]. As a consequence we shall bound from above (§5) the 
spacing between invariant tori, i.e. the gaps sizes, and apply the results to show that 
the homoclinic angles are (generically) very large compared to the (dimensionless) 
spacing: this will imply (§6) the existence of long heteroclinic chains for a system 
consisting of a pendulum interacting with two rotators thereby deducing (by using 
§8 of [CG]) the existence of drift and diffusion in the system when the perturbation 
parameter e is smaller than a power of n, rather than exponentially small in r]~^. 
Finally we show (§7) that the perturbation can have a large component which, 
provided it depends only on the fast rotator angle (or more generally provided it is 
"monochromatic", see §7), does not alter some of the results. 

A problem with three time scales naturally arises in celestial mechanics: in the 
theory of the precession motion of a satellite axis due to its non sphericity. Under 
certain approximations, see [CG], it is described by a Hamiltonian of the form (1.1) 
plus a large "monochromatic" term. In [CG], §12, it was shown that Arnol'd diffusion 
occurs, but a preceding error on the lower bound of a homoclinic splitting infirms 
the conclusions. In [GGM1] the correct lower bound is given and combining [GGM1] 
with the new results of §7 one sees that the analysis in [CG] applies, provided that the 
monochromatic term is absent (this means that the normal form method described in 
§8 of [CG] applies to the heteroclinic chains that, below, we show to exist). Finally, 
in the last section, we outline a possible strategy on how to approach and complete 
the analysis of the full problem of §12 of [CG]. 

The non standard part of this work relies strongly on the explicit formula for the 
homoclinic splitting derived in [GGM1] and called the ll large angles theorem" . 

§4. Averaging. 

We shall first perform a canonical transformation, that we shall denote to 
"reduce" the perturbation to order 0(e^), with Af an integer suitably large, but 
such that no "small divisors" appear in the definition of $/v". This is easily achieved 
by requiring N < (n~ 1 N~ 1 , with £ = |u>2|(2(D(,) _1 , (note that the existence and the 
relative sizes of three time scales is essential here) . Subsequently we shall apply the 
results of [CG] to the Hamiltonian expressed in the new variables; in this way we 
shall be able to prove a more refined bound on the size of the gaps between invariant 
tori, (and show that it is smaller than the homoclinic splitting without the need to 
require e to be exponentially small in n~^). 

Fixed an integer number J\f, we look for a function $>s/(r]%A[, a,p',q) generating 
a canonical transformation which casts the Hamiltonian (1.1) into the form (3.3) up 
to order A/". This means that, in the new variables, the Hamiltonian will have the 
form (3.3), but with a different function G (that we shall denote GV) and an extra 
additive term (the new perturbation) which will not be of order 0(e) but of order 
0(e M ). 

The functions <I>w, G_\r will be polynomials in e: 

G N =eGV+s 2 GM + ... + ^GW (41) 

We proceed to determine recursively the coefficients for N < 7}'n~ 1 N~ 1 in order to 
prove the following averaging theorem. 
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Theorem I. Fixed Af < Cv' 1 ^' 1 , with ( = |w 2 |(2J) h )- 1 ; the functions $0'),G/0) 
exist for all 1 < j < Af and are holomorphic in the "decreasing" domains Tf> Pj ^ jyKj 
with £j = £ — jd, Pj = poe^i 5 , Kj — Ko e ~ : ' 5 f or 5 > such that MS < i£ ; an d they 
verify the bounds: 

|$0)| < DB j -\j - 1)!, \GW\ < EoB'-^j - 1)!, j < Af (4.2) 

where D = bEoT^Do and B = b (EqTq 1 X^ 1 tl^ 1 ) 2 D 2 Af 2 , if the constant b is large 
enough, and D = ^^^{Af^ 1 ) 2 log^C^ 1 ). 
Furthermore: 

• Map (3.2) with $ = $_a/ generates a canonical transformation from ^ 
f (0o e-«o/2 i£o /2, Koe -«o/2 to the domain X> Po , ?0 , Ko if: 

6' sma/Z enough. 

• ,4nG? in t/ie new variables (A', a',p',q') the Hamiltonian (1.1) becomes: 

h(^A[) + v -iw 2 A' 2 + G( P 'q', rM) + Gjtf(pV> ?M)+ 

/ e \^ i ( 4 - 4 ) 

+ (^) /v, e (»7*4, a'.pV) 

with |e- 1 G^-|,|/v, e | < 2S . 

The proof follows immediately (for instance) from the arguments used to prove 
Nckhoroshev's theorem in [G1,BG] and is given in Appendix Al, for completeness, 
extending similar considerations in [CG] . 

Remarks. (1) An interesting consequence is obtained by choosing Af = jn^i with 
7 > 0. In this case under a condition |e| < e = 0(r)2 / log 2 r/ -1 ), see (4.3), the 
Hamiltonian (1.1) can be put in the form (4.4). 

(2) Another interesting consequence (of the proof of the above theorem given in 
Appendix Al) concerns the case with Af fixed independently on n and a perturba- 
tion / whose Fourier transform with respect to a does not vanish only for wave 
vectors y_ multiples of a given wave vector v_ {"monochromatic" , or "integrable" , 
perturbation). In the latter case the condition is |e| < 0(1) for all Af, prefixed; see 
the corollary discussed in the last paragraph of Appendix Al, (although it would be 
easy to prove directly the simple statement made there). 

(3) As a quantitative estimate of the parameters of the new Hamiltonian (4.4) it 
follows from the proof in Appendix Al that the parameters Xq,T , E , M can all 
be taken the same as those of the initial Hamiltonian within a factor 2: also n is the 
same while, of course, e becomes (e(2£ ) _1 ) A ^- 

§5. Abundance of tori. Interpolation. 

Given the Hamiltonian (4.4), define the free spectrum So as the collection of the 
rotation vectors 

uj = u(x',r)i A') d = f (nid'hiniA'J+nid'Gix'^iA'J^-iuji) (5.1) 

as A[ varies in [— Rn~i , Rn~i] and x' = and d' denotes differentiation with 
respect to the argument a' = 772 A[. More generally we define S to be the collection 
of vectors (5.1) as A[ e [-Rn~i , Rn~^] and \x'\ < K e^^'/ 2 . 
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Let us consider the values A[ such that the rotation vector u G So in (5.1) satisfies 
at x' = the diophantine condition: 

+ 77^2^2 > C(r/)|^r r for all ^ v e Z 2 (5.2) 

where //'Wi = t]?loi(x', rj^A^) is the first component of o> in (5.1), and r > 1, 
C(ry) > are Diophantine constants: ui has size of order 1 (see §1). 
Such A[ occupy a set in [— Rr]~2 , Rr]~ 3 ] whose complement has measure bounded 
above by: 

b"M E Q \^C{ri) (5.3) 

with a suitably large b" , because the measure of the set of cji's in the interval 
of variability of Wi ([cD a , u>b] C [0, SoAq 1 ], see §1) verifying (5.2) has complement 
bounded by b"E \g 1 r]^C(rj)), for some constant 6", and the derivative of the map 
from action to frequency is T]d 2 h(pq,r]i A{), bounded below by nM^ 1 (see (2.2)). 
Hence any interval in [— Rn~ 2 ; Rrj~ 2 ] of size A(iy) > b" M E \Q 1 C(r]) will contain 

points verifying (5.2). 

_ 1 

If C(rj) = Qe~ sn 5 with s, ft > (the factor f2 is here just to fix the dimensions 
and it could be any constant with the right dimensions, e. g. LO2 or Eq/Xq or other), 
then it follows from the combination of theorem 1 above and lemma 1' of [CG] that 
for each lo € Sq verifying (5.2) the Hamiltonian system (3.3), hence (4.4), will have 
an invariant torus such that the rotation vector lo ' assumes the value lo , for 77 small 
enough and |e| < eq. 

In fact the "effective coupling" , by applying theorem 1 with Af = 777^2 ; is of order 
_ 1 

O((e/2eo) ln 5 ) and, for 7 large, this is so small that even if multiplied by {C(n)/n]~ q 
for q > arbitrarily fixed gives a still very small result when r\ gets small enough. 

The stability of the tori and their stable and unstable manifolds or of the motions 
in their vicinity depends on the size of the latter product (for a suitably large q, see 
[CG], cq. (5.76), where q — 6 is an estimate) which is small if 7 > gs(log2) _1 . 

Note that if we simply applied lemma 1' of [CG] to the Hamiltonian (1.1) we would 
have found, as a bound on the convergence radius for the parametric equations of the 
tori and of some of their nearby motions, a quantity of order O(so) times \C{rf)/VL] q , 

which would have been exponentially small in 77 2 if C(rj) — fie" 1 "' 2 with ft, s > 0. 
Another way of interpreting the analysis (see [CG]) is by saying that there is 
a change of coordinates (A, a,p,q) < — ► (A', a',p',q'), defined in the vicinity of 
[Refill < r/~2_R and \q\ small, which is of class C°° (in (A', a',p',q') and e) 
which casts Hamiltonian (1.1) or (4.4) in the form (3.3) on the set of (A 1 , a',p', q')'s 
for which the vector lo_ ' = L0_'(p'q' , fp^'i) defined in (3.4) verifies (3.5), i.e. it is equal 
to some w G (So verifying (5.2). 

Furthermore, fixed A' and p',q' in such a set, the coordinate transformation (of 

_ 1 1 

the remaining two coordinates a) is analytic in e for |e| < eoO({C(i])/£l} q ' r 1,5 ). 

And the identity between (1.1) and (3.3) with (3.4) holds, on the set on which the 
functions agree, also between their first derivatives with respect to the (A 1 , a',p',q') 
variables. So that the sets with A' and p'q' = x' fixed and lo' verifying (3.5) are 
invariant and the motion on them is very simple, and described after (3.3). 

We stress again that the latter statement is a slight generalization of the quoted 
results in [CG]: namely the latter really covers the statements in the preceding 
paragraph for vectors lo_ '{x 1 , rp A^) which verify (3.5) and have the special form 
w'(x',?7^Ai) = (1 + 7(2/, rjiA[)) w'(0,77^Ai) for some small 7(2;', r)i A[). This 
restricts us, de facto, to considering cases less general than those that verify (3.5). 
However, as mentioned above, if one goes through the proof of lemmata 1,1', 2 of 
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[CG] one realizes that the proof covers, unchanged, the more general case we quote 
here. 

As a consequence the following high density theorem of invariant tori also holds 
(corollary of theorem 1 and of the above analysis). 

Theorem 2. Fixed 0,s > let C(r?) = fle- sr >~* ; and let JV = 777-3, 7 > 0. For 
r) small and \s\ < So, where So is given by (4-3), the Hamiltonian (1.1) has, for 
all 77 small enough, a family of hyperbolic invariant tori, with stable and unstable 
manifolds of dimension 3, whose rotation vectors lu' = (772^, t?"^) G So ver- 
ify (3.5) and therefore have first components which fill the interval [w ,u>(J within 
_ 1 

0(£le~ srl 2 ) for all 77 small enough, provided e < Eq. 

More precisely the parametric equations of such tori can be written, for e small 
enough and u' e So verifying (3.5), as: 

A=H^,{a>) I =!„,(<*') 

QL =a'+h„,{a') ¥> =*«'(«') 

with H_ u >, h u i,Iu'^ u'> at fixed lu', analytic in e and in a' E T 2 and divisible 

by e, so that such functions can be found by perturbation expansions. The small- 

_ 1 1 

ness condition on e is \e\ < 0(eo[C(r])] ql T ' 2 ") (a very weak condition, as 7 can be 
arbitrarily prefixed) and the motion of the data (5-4) is simply a' — > a' + Lu't. 

It is interesting to write the complete condition that e is small enough (including 
the constants) and a more quantitative expression for the filling of the action axis. 
Since the above theorem relies on [CG] one has to find bounds for several quantities 
associated with (4.4). If Mq denotes a lower estimate for the minimum value of 
the second derivative of (1.1) evaluated at e = with respect to a — 772 A\ {i.e. if 
M is defined as in (2.2)) one finds, applying the condition (5.90) of [CG], for some 
B,q>0: 



kl < ^0 ^gL €o C(r?) ) q,U = const . eo [C^r- 1 ^ = 0(77!+) (5.5) 

having unified various constants to simplify the expression (and after translation 
to the present symbols and conventions) and denoting |+ a prefixed number big- 
ger than |. And under the same conditions the filling of phase space is within 
0(Moi?oA 1 C(77)). 

A more detailed description of the above analysis is the following strengthening of 
the above theorem. 

Theorem 3. Under the same hypotheses of theorem 2 above there exists a C°° 
function T~to(A ,x'), x' — p'q', defined in 2?p /2,£ /2,K /2 an< ^ a C°° canonical trans- 
formation, defined on the latter domain, and having the form: 

A=A'+ H(A',a',p',q% I =2(A', a',p',q') 
a=a'+h(A',a',p',q'), ip = $>(A', a',p',q') 

with H_, h,!,^ of class C°° divisible by e and analytic in e, a' at fixed A',x' = 
p'q 1 such that if lu' = d A 'Ho(A' , x') verifies (3.5) then in the new coordinates 
the Hamiltonian coincides, together with its first order derivatives, with l~to(A', x'). 
And Ho is close within 0(e) to h(q^ A[) + 77" "2 uj^A^ + G(p'q' ,772 A[). Hence the 
motion of the data (A', a',p',q') for which u' verifies (3.5) are A' — > A', a' — > 
a' + u't, p' -^p'e-9'\ q' -> q'e+s't with g' = d x ,Uo{ A ', x') ifx'=p'q'. 
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Remarks. (1) The value C(n) = Qe sn 5 is particularly interesting in view of the 
results in [GGM1], where the Hamiltonian 

i A 2 i I 2 

H = r)*A x + n-± + t]-^A 2 + — + 5 2 (cos tp - 1) + e/O, ai, a 2 ) (5.7) 

is studied, with / an even trigonometric polynomial of degree N, Nq in a , tp respec- 
tively and the homoclinic splitting of the invariant tori whiskers is found to be of 
size of order C(n)\ s=7r / 2g - Note that (5.7) with (7, tp) near (0, 0) is a particular case 
of (1.1). One just uses "Jacobi's coordinates" (p,q) instead of (I,<p), see Appendix 
A9of[CG]. 

In [GGM1] the slow frequency is considered to be i] a , with a > 0, so that a = 1/2 is 
only a particular case. Of course in what follows we could also consider an arbitrary 
value a > 0, and essentially nothing would change, but we have preferred to confine 
ourselves to the case a = 1/2 for definiteness. 

(2) In general the functions H , h in theorem 3 do not have zero average for p' = 
q' = 0. The even symmetry of / and of the pendulum Hamiltonian in (5.7) imply 
that the variables A' have a simple physical interpretation for the points on the 
invariant tori: they are just the time averages of the unperturbed actions. In fact 
this " absence of torsion" is a very remarkable and useful property of the system 
(5.7) (Thirring model) which was pointed out in [G3,G4] (and called the property of 
the tori of being twistless). It is due to the special symmetries of the system (5.7) 
and to the separation of the energy into a quadratic part involving actions only and 
an angular part involving only the angles. 

(3) A further consequence of the symmetries of (5.7) is that each invariant torus 
described by theorem 2 has stable and unstable manifolds that intersect along a 
trajectory, a homoclinic trajectory, that when tp = n has angular coordinates a = . 

The difference at a,tp, that we denote Q ( a , tp; Tp A' x ), between the A -coordinates 
of the stable and unstable manifolds of the torus T(A' 1 ) with average action A[, 
vanishes at tp = tt, a — 0, sec [CG], §9. Hence we say that a = 0,tp = tt is 

a homoclinic point for T(A[). Setting Q(a) d = Q ( a , tt; n^A^) the matrix = 
d ai Qj(a)\ a = o is called the intersection matrix and in general has a determinant 
det D that does not vanish: its value is called the homoclinic splitting, sec [GGM1]. 

(4) The same result as theorem 2 can be obtained, for the Hamiltonian (5.7), also 
from the Lindstcdt's series expansion for the stable and unstable manifolds of the 
invariant tori, by using the methods developed in [E,G3,Gcl,Ge2] and an (improved) 
version of Sicgcl-Bryuno's lemma; but this is a completely different approach to the 
problem and we expose it, in detail, in a forthcoming paper, [GGM2]. The new 
approach yields, as expected, a better bound on the convergence radius 6q: namely 
eo = 0(n 2 ) compared to the ~ 0(?7 4 - 5+ ) that would be obtained by applying theorem 
1 to (5.7) (see the remark 1 after theorem 1 and the comments after (5.5)). 

(5) Theorem 3 means that there is a canonical system of coordinates in which the 
Hamiltonian takes the normal form Ho given by an expression like (3.3) on a set 
with very small complement, in general not empty and open (we stress that this is a 
statement different from the one that would hold if the Hamiltonian could have been 
really put in the form (3.3), i.e. in the whole vicinity of the real domain of (1.1)). 
The reason Tlo is called a "normal form" is that the data with uj' verifying (3.5) 
evolve in a very simple fashion i.e. as A = const., a' — > a' + w'f, p' — > p'e~ 9 1 
and q' -> q 1 ^' 1 (see (3.3) and (3.4)). 

(6) Another consequence of theorem 3 is that fixed A' so that oj' = d A>T~io(A' , 0) 
verifies (3.5) it is possible to fix Ai(x') so that d A''rto{A l '(x'),x') = lo_' for, say, 
\x'\ < ko/4. And, at the same time, the energy of the motions that start in 
£po/2,£o/2,K /2 witn A' = A'(x') and p',q' with p'q' = x' and with a' G T 2 
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arbitrary is x'-independcnt. Here the linearity in A' 2 is used to adjust the energy 
once A[(x') is determined: so that the apparently useless role of the variable A' 2 , 
"reservoir energy" for keeping the second rotator in a constant speed rotation, can 
be well appreciated here. This is a key property for the theory of diffusion in phase 
space developed in [CG]. 

§6. Existence of heteroclinic chains. 

The novelty of the results in §5 is the possibility of using the "large angles theorem" 
of [GGM1] in the same way as the erroneous result (with the same name) in §10 of 
[CG] was used to discuss heteroclinic chains. 

In this section we restrict considerations to the system described by the Hamiltonian 

(5.7). Given an invariant torus with average action A such that the corresponding 

_ i 

rotation vector uj_ verifies (5.2) with C(n) = Qe~ sn 5 , O, s > 0, one has that the 
torus and its stable and unstable manifolds (whiskers) are analytic for |e| < 2e , for 
a suitable £q. The discussion in §5 gives £o = 0(r]^ + ) (by taking J\f = 777^2 with 7 
large enough, see (5.5)), by the argument leading to theorem 2; as mentioned in §5, 
remark 4, a better value £0 = 0(rj 2 ) can be obtained with a different method. 
Then we want to infer that the homoclinic splitting, at the homoclinic point 
(<p, 01,02) — (tt, 0,0) (see §5, remark 3), of such a torus is generically given by 
the first order (Melnikov's integral), i.e. it is of the form: 

CTeV'e - * 1 " 819 " 1 ""* (6.1) 

with a a non vanishing constant and b a positive constant linearly depending on the 
degree Nq (in ip) of the perturbation /. The finiteness of No is thus a key assumption 
for the validity of the following analysis. 

The theorems of [GGM1] do not apply directly because in the latter paper C{r}) was 
taken proportional to n d for some d > 0. On the other hand, by the considerations 
in §5, the results of [GGM1] can be immediately extended to cover also this case, 
i.e. to conclude that the splitting is given asymptotically by Melnikov's integral, as 
we shall see below. 

By theorems l-j-3 above, the invariant tori and their stable and unstable manifolds 
are analytic in e for |e| < 2e , with e = 0(n c ), where d > 9/2 or c' — 2 (according 
to which bound for e is chosen, see the remark 4 in §5). This implies that, to any 
order h, the parametric equations of the manifolds (i.e. h,H_&s functions of a,e 
at x 1 = p'q' = 0, A[ fixed) evaluated at ip = ir are analytic in e and their Taylor 
series coefficients can be bounded to order h by: 

D 2 B , i- 1 (e2- 1 e^ 1 ) h (6.2) 

for some positive constants D 2 , B 2 . 

Moreover we know, from the analysis in [G3] and [GGM1], that, to any order 
h < Af = 77?~ 5 , the Taylor coefficients for the splitting are bounded, to order h, by 

D s h^B^(sr,-P') h e-^\ 9 ~^~ h (6.3) 

with = 2(Nq + 1) + 2, for some positive constants D3, B$ and with No equal to the 
degree in p of the perturbation / (hence 7V < N: the convenience in distinguishing 
between N and TVo is that it is possible to extend the present work to cover cases 
in which / is analytic in a, i.e. N = 00, provided one keeps it polynomial in the 
angles ip). 

No factor C(t]) appears in (6.3) as no resonances occur in the bounds of the divi- 
sors that appear in the perturbation expansions (of [G3,GGM1]) contributions to the 
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splitting at orders < Af, by the very choice of Af. In fact as long as only orders < Af 
are considered one has | w ' ■ v_ \ > 772 as \v \ < N-fn~z = NAf. 
Then if e < 0(n c ), c > max{c',/3' + 2} = /?' + 2, the contribution to the split- 
ting of the orders > Af can be estimated from (6.2) by O{(e2~ 1 e~ 1 ) A/ '): which is 
much smaller than the generic homoclinic splitting, (6.1), evaluated from Mclnikov's 
integral. 

The contribution arising from the orders 3 < h < Af is bounded, via (6.3), by 

Oiierj-Pfe-^^ 1 ^^), with /? = /?' + 2, see [G3,GGM1], also much smaller than 
(6.1) if e is smaller than a high enough power of 77 (i.e. without a condition of 
exponential smallness of e). 

Hence we conclude that under a condition like |e| < 0(i] c ), with 77 small enough, 

tori with average action A , arbitrarily prefixed energy and verifying (5.2) with 

_ 1 

C(r}) = f2e~ sr ' 2 , 0, s preassigned) exist and their homoclinic splitting has size 
_ 1 _i 

0(e 2 r/^ b e^i^ 2 ^ 9 v 2 ) generically in perturbations / which are trigonometric poly- 
nomials in ip (and the constant b depends on the degree Nq in ip of /). 
For instance the simplest case (5.7) with / = cos(ai + ip) + cos(a 2 + <p) and g = 1 
gives a leading splitting, see [GGM1] eq. (6.14) and (7.19), as 77 — ► 0: 

32irr ) -h 2 e-% n ~ i (6.4) 

We also conclude the (generic) existence, for s < 0(n c ) and 77 small enough, of 
chains of isoenergetic invariant tori ("heteroclinic chains") with average actions A'^ 
such that: 

(i) w '(0,77-^1) verify (5.2) with C(t?) = Oe" 8 ''"' with s,fl that can be fixed a 

a - 1 

priori, and | ^4^- — Aj+il < Oir\~ c e~ sv 5 ) for a suitable c"; 

(ii) the homoclinic splitting has size 0(e 2 77~ fc e _ ^'" 2 ' s ^ 5 ) (i.e. far larger than 
C(?7) if s is fixed a priori to be s > fg _1 ); 

(iii) the unstable manifold of the torus with average action A j intersects the stable 
manifold of the torus with average Aj +1 in a heteroclinic trajectory with splitting at 

p = 7r of order 0(e 2 rj^ f) e^i I" 2 ' 9 ^ this follows, as usual, from an application 
of the implicit functions theorem because the equation to be solved for finding a 
heteroclinic point is an implicit equation with Jacobian determinant, at the trivial 
solution a = (corresponding to the homoclinic point at <p = n of A'j), given by 
the splitting detD, see also [CG.GGMl]; 

(iv) the "genericity" is a very explicit condition because it simply requires the lowest 
order value of the splitting to be non vanishing: this holds generically in / (picked 
up inside the class of functions we are considering); 

(v) furthermore, since the set of points verifying (5.2) is abundant in the sense of the 
above theorems 2 and 3, we also obtain (at no extra cost, i.e. by "abstract reasoning" 
on measure theory) that the values of A'^ x are density points for the set of values 
A'ji £ [— r]~iR, T]~iR] whose corresponding rotation vectors verify (5.2) with the 
chosen value of C(?7); and we can suppose that the values of A'^ at the extremes of 

the chain are close to the extremes of the interval [—i]~iR, n~i R] within e (or even 
closer) . 

We can summarize the above discussion in the following result. 

Theorem 4. Given the system described by the Hamiltonian (5.7), there is a se- 
quence of invariant tori with average actions A^, . . . , A^, with \A' 1 — A'j^\ ~ 

2Rr]-h (within 0(e)) and \ A' i - A' i+1 \ < Re-^^ for i = 1, . . . ,7V - 1 and s > 0, 
such that the unstable manifold of each invariant torus intersects the stable manifold 
of the torus following it along the chain. At the heteroclinic point at ip — n the 
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splitting between the two manifolds is given by a' e 2 n b e Sl^g v 5 ; where a',b 
are constants depending on the perturbation f, with a' generically not vanishing. 

This extends the theorem of existence of heteroclinic chains discussed in [GGM1], 
§8, to the anisochronous case. 

§7. Fast averaging theorem and abundance. 

A second remarkable application is a fast averaging theorem which follows from the 
proof of theorem 1 and concerns a system with Hamiltonian: 

i I 2 

H = h(n^A 1 )+n-^ 2 A 2 + — + J a g 2 (cos(p-l)+eifi(a,(p)+e 2 f2(a,tp) (7.1) 

ZJq 

where I,ip are conjugate variables, Jo, 9 are strictly positive constants and £i,£2 
are perturbation parameters. The function f\ is monochromatic (i.e. with Fourier 
transform not vanishing only on modes different from zero and multiples of a fixed 
u_ = [yx^v-i) with v 2 ^ 0), while f 2 is only required to be a trigonometric polyno- 
mial. 

On the basis of the results of the previous sections we can expect that in order to 
have "long" heteroclinic chains of invariant tori we need to require e\ and e 2 to be 
both small (of the order of e of theorems 1-^4). 

But this is not the case because of the monochromatic nature of /i and we shall 
show that the results of theorem 1-^3 hold essentially under the assumptions that 
e 2 is small of order e while E\ can be quite large, up to almost 77^2; in particular 
the value z\ = 1 is amply allowed. In a sense also theorem 4 could hold; however we 
can only show (see below) that verification of its validity can be reduced, at least in 
principle, to a finite computation if 77, e 2 and E\ — 1 are fixed with e 2 < Eq. 

We consider first the case e 2 — 0. In this case, by performing the canonical trans- 
formation with generating function S = A' ■ a + I'tp — £\rj^ oj 2 d 2 fo( ®_ , </?), one 
sees that the size of the perturbation can be reduced to 0(n 2 ei) under the condition 
that r\^E\ is small enough (so that the canonical map generated by S makes sense 
in a domain slightly smaller than the domain of h). 

This is remarkable because the change of coordinates is globally defined and it is 
not restricted to the vicinity of / — ip — 0, and it leaves the perturbation still a 
monochromatic trigonometric polynomial with zero average over a . This cannot be 
pursued to higher orders without losing the trigonometric polynomial nature of the 
perturbation. 

Then we can then use Jacobi's coordinates near I — <p = 0, see [G2], to put the 
Hamiltonian in the form (1.1) with sf of order 0(n^E\). 

It follows from the proof of theorem 1 (see the remark 2 in §4) that if the pertur- 
bation is monochromatic the condition for being able to cast the Hamiltonian in the 
form (4.4), with a prefixed and 77 independent Af, say Wo, is that \ei\ < 0(1), so 
that we see that (7.1) can be cast in the form (4.4) with a prefixed Afo for \ei\r)? 
small enough. 

We consider now the case s 2 ^ and we imagine performing the change of variables 
that puts the part of (7.1) not containing the last term into the form (4.4) for some 
Wo. 

This casts the full Hamiltonian (7.1) into the form (1.1) except for the fact that the 
perturbation is no longer a trigonometric polynomial in the a variables: it is just 
an analytic function in a domain of size £0 > in the complex planes for the a 
variables. 

We truncate its Fourier expansion in the a at \ v\ < N, where N is left as a 
parameter. The remainder can be bounded in its domain by ae~^« N for some a > 
of order Odeirp | + |e 2 |)- 
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We arc, therefore, in a position to apply again theorem 1 and to cast the full 
(7.1) into the form (4.4). We make the following choices: N = r\~ 1 -\J (logry -1 ) -1 , 
N = \f)~% yf\ogrf~^, (assuming r\ < 5, which is not restrictive as we are interested 
in properties holding as r\ — > 0). 

Combined again with theorem 1 above, this means that we can put the Hamiltonian 
in the form (4.4) with the last term replaced by an "effective interaction" of order: 

/ (( £l ^)^+ £2 )(logr ? - 1 )-| y^Qo g .- 1 )- 1/2 | ac _^ v -i aoKV ^ ^ e _ zv -i 

(7-2) 

provided £2 — 0(r) c ), with c large enough and A/"o so prefixed that 7Vo/2 > c; a,£o 
are constants that (if wished) can be explicitly estimated from the above argument. 
The inequality holds for any prefixed z > provided -q is correspondingly suitably 
small. 

This means that, in some sense, the size of the monochromatic term in (7.1) does 

not matter too much and e\ = 1 (or even almost £1 ~ 77^2) is a sufficient condition, 

together with |£2 < £0 and rj small enough, to guarantee the existence of an invariant 

torus which is run quasi-periodically with a rotation uj' verifying a (very weak) 

_ 1 

Diophantine condition {C(rf) — Q,e~ sr * 5 ) and whose equations can be constructed 
by a convergent perturbation series. 

A further consequence is that the Hamiltonian (7.1) has many invariant tori whose 
average actions fill the action space as described in theorem 2 above, so that we can 
apply the results of [GGM1]. We summarize the above discussion in the following 
averaging theorem. 

Theorem 5. Given s > consider the Hamiltonian (7.1) with fj verifying the 
assumptions following (7.1). If \s\\ < r\~ d , \e 2 \ < V +d w ^ prefixed < d < \ and 

d' > § there is a family of invariant tori T(A') with average actions A' filling action 
_ 1 

space within 0(e~ sv 5 ). The parametric equations of such tori can be computed by 
perturbation theory (i.e. by convergent power series in £1,62 J o-rid they are, together 
with the homoclinic splitting, analytic functions of £\,Ei provided r\ is small enough 
( depending on the size of s,d, d! ). 

Note that the "torsion free" property allows us to fix the rotation vector of the 
motions on an invariant torus by fixing the average value A\ of the action, no matter 
which is the perturbation size e: in fact in the model (7.1) that we are considering 
it is lo_' — (772 + r\A\,r\~?-i02) which is £-independent (remarkably enough). 

Apart from the intrinsic interest of theorem 3 it has some relevance for the theory 
of drift and diffusion in a priori stable systems, like the precession problem whose 
analysis was attempted in [CG], §12. We discuss this aspect in the next section. 

§8. Drift and diffusion under large monochromatic forcing and small quasi 
periodic perturbation, or in a priori stable systems 

The interest of the above theorem 5 is, again, that it might be used, by repeating the 
argument in §6 to show the existence of heteroclinic chains with the same properties 
of the ones discussed in that comment even if E\ = 0(1). 

This time there is, however, an extra difficulty if Si is large, i.e. really of order 1: 
the argument of §6 provides us with a tool to bound the remainders, but not to find 
out which is the dominant term among those of "small" order (i.e. smaller than J\f 
in §6). 

The reason is simply that the contributions of higher order in £1 not only are not 
smaller than those of order 1 but in fact are checked to be a priori larger. Meaning 
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that one can find some contributions to them coming from suitable graphs, in the 
sense of the representation in [GGM1], and larger than the lowest order contribution: 
but in a subject where cancellations are quite common this does not seem to mean 
much. 

However the results in [GGM1] and theorem 4 above show that the problem of 
proving the existence of a heteroclinic chain joining the extremes of the action A[ 
interval [—Rr]~^,Rr]~^] can in principle be solved by a finite calculation, if e\, £2 
and 77 are fixed. The program is illustrated in the following. 

In fact such chains can be proved to exist and to consist of finitely many elements 
for £i| small enough (by theorem 3) and the above discussion shows that the series 
in £i,£2 for the splitting is analytic, and the term of order hi in e\ and hi in £2 is 
bounded by 

B hl+h2 (s 1 r)i) hl {e 2 ri~ c ) h2 (8.1) 
where B is a suitable constant. 

Assume that we can show that the sum of the first hi+hz < 777^2 terms in the series 

- 1 — - 

for the splitting is bounded from above and below by a quantity 0(e~ 2 l w 2|g v 2 ); 
then the remainder is surely smaller than this quantity. Hence we can approximate 
the heteroclinic points for an allowed value of £\: keeping fixed the average action 
A[ we can increase it following the evolution of the heteroclinic point. This is a 
finite calculation, in principle, because we start from E\ > and have convergent 
expansions up to £1 = 1 (and beyond up to almost r)~? ) thanks to the above fast 
averaging observations, valid close enough to the tori (where the changes of variables 
discussed to reach (7.2) make sense). Thus with a finite calculation one can get 
controlled approximations of the values of the heteroclinic angles (at fixed £1, £2, 77). 

What is not a priori guaranteed is that such angles, which start positive and "large" 
compared to the separation between the tori, do not become small while £1 grows. 

However there is no reason for this to happen: one expects the angles to become 
wider, not smaller, or to vanish only for finitely many values of £1, at least generically. 

Hence a computer assisted analysis of the problem is likely to succeed and it might 
even be used to reach values of rj and £i,£2 where the angles are large and the 
spacing between the tori can be widened (thus diminishing the number of invariant 
tori of the chains) . As discussed in [CG] the drift along a heteroclinic chain with non 
zero splitting is always present and one would obtain models with Arnol'd diffusion 
in systems like (7.1) with £1 = 1. This would correct even the applications of the 
incorrect version of the large angle theorem of [CG] (see §12 of [CG]) and make it 
rest entirely on the corrected version in [GGM1]. 

In fact a problem like (7.1) with £1 large (£1 = 0(1)) arises in the theory of 
Arnol'd diffusion in a priori stable systems, like the precession problem treated in 
§12 of [CG]. From the analysis in [CG] it seems that such large (but monochromatic) 
coupling Hamiltonian systems might arise often in the reduction of a motion near a 
multiple resonance of a a priori stable system. 

Remark. R is interesting to note that the above program with h(a) = W\a, i.e. (7.1) 
in the isochronous case, can be carried out without need of a computer as remarked 
in [GGM1J (the diffusion then follows via the results in [G2]) at least for all E\ < 
0(n~ 2 ) except possibly a finite number of them. 

Acknowledgments. This work is part of the research program of the European 
Network # ERBCHRXCT94 04 60 on: "Stability and Universality in Classical Me- 
chanics ", 
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Appendix Al. Proof of theorem 1. 

Suppose that for h < k — 1 one has: 

< DB h ~ 1 (h - 1)!, \\G (h) \\ < EoB' 1 - 1 ^- 1)! (Al.l) 

where D is such that the bound holds for h = 1. Here || • || denotes the maximum 
in 'D ph ^ h ,K h Wl th £/, = £o — ^>P?i = Po e ~ Sh , K h = K a e~ hS . Furthermore the Fourier 
transform of $W with respect to the a -variables has only modes i/ with | £ | < AT/i. 
Therefore for h < (r/~ 1 N~ 1 , with ( = \u 2 \{2u) b )^ 1 , we can bound the various terms 
in the difference between the two sides of (3.1) as follows. 

We expand in Taylor's series the functions in (3.1) with respect to the point ob- 
tained by setting <I> = 0. We define u d — (rj? dh(r]i A[) + r/idG(p'q,r]^A' 1 ),ri~?uj2), 
gip'q^^A^) = d x G{x,r]^A' 1 )\ x=plq , d = qd q - p'd p > and A d = u -da + 90- Then 
the linear term in $( fe ) is simply A& k \ 

The nonlinear terms arising from h(T}^(A[ + <9 ai $)) contribute to $( fe ) an amount 

(k) 

X[ such that: 

Axf^f^^M" E n^^)* (^1-2) 

a>2 k 1 + ... + k a = k j-\ 

kj >l 

where the * means that from the Taylor expansion in p' , q and Fourier expansion in 
a of the term in parentheses one has to subtract the terms which have Fourier 
mode and that have an equal power of p' and q. In fact the latter terms make up the 
contribution to of the nonlinear terms arising from h(r]2 (A[ + d ai $)). After 

(k) 

this subtraction the operator A is invertible, i.e. X[ is well defined and the action 
of A -1 can be bounded by 1 C(r], S) because as long as \ u\ < (r/^ 1 one has: 

+ ki^i+i^i+h+h) 

E ,, | rr. i -in < r o ^(V, 5) (A1.3) 

u, n ,m (H 7 ? 2 + I^W 2 +\m-n\)T 

where the + means that the sum runs over the u_,m,n for which the divisor does not 
vanish and, furthermore, \ u\ < (r/^ 1 , so that no "resonances" occur. The quantity 
C(rj, S) can be taken 

C(rj,S) =c(ji*5- 3 +r)~ i S~ 2 iog^ 1 > 1 ( A1A ) 

where c is a suitable constant (recalling also that e, rj < 1, see §2). 
Making use of the analyticity assumption, which implies that \h^\ < E p^ a in 

l^l<^E(Tr) a ^ E ftte-D" (^-5) 

a>2 ^° fci + ... + fc a =fc j=l 

kj >l 

where the factor 2 above comes from the estimates of the Fourier transform of the 
factors d ai $( k ^ in (A1.2) inside the domains D Phj+1/2 ^ kj+1/2 ^ k]+1/2 . 

We shall repeatedly use below the inequality: %+.. .+k a =k rii=i(^j — !)■ — — 1)' 

kj >i 

(to prove it, simply bound the product of the factorials by the (k — a)\, and the 
number of addends by (k — a + l) a_1 /a!). Then, from (A1.5), 

^<- E ^-^& «%£<\ <-> 
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Hence in D Pk . Sk . Kk : 

The higher order terms contribution from the third term in (3.1) is: 

n a 

Ax 2 k) = ( e G(a,n) e ^"n^^n^* ) (j4l8) 

»•<• + + j = l j=l 

a + n>2 f ll + ... + h„ = fc 

with hj,ki>l. So that in D Pk ^ ktKk : 

|x? , |s ^ 5 ^^-„-^ (t _ 1) , £ 

a + „>2 

provided: 

_2D_ 1 4 W) 4P 1 

«5A i? 2' r A5W < 5 1 j 

and note that (ALIO) implies the conditions in (A1.7). 

Likewise the fourth term in (3.1) yields a higher order contribution bounded 
as: 



LA1.9) 



^ ^V^J" (ALII) 

v y L A B<5 A B6 2 

so that one has in D Pki £ ktKk : 



The fourth higher order contribution comes from the third term in the r.h.s. of 
(3.1) and is bounded in the same way as leading to the same bound under the 

same conditions. The fifth higher contribution comes from the last term in the r.h.s. 
of (3.1), which, up to the term & k \ is bounded in D Ph £ htKk as: 



|*f I < E *o n B k -2»D»±(k 1)! < 



n>l 



< 2 EoC{71 ' 6) 2D B k -Hk 1)! < l -DB k ~\k 1)! 
Lo XqBS 5 



(A1.13) 



provided: 



2D < 1 2 E C( V ,S) 2 < 1 (i4LM) 



A B(5 2 L A .B<5 5 



It remains to estimate itself. Noting that it is the collection of the terms 
that are subtracted by the above * operations it is clear that it is bounded by 
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DB k ^ 1 (k — l)\T()C(ri, S)^ 1 because one does not have to invert A to find it from 
the bounds on the higher order terms of (3.1), and one thus saves a division by 
r C(?7, (5) _1 . Hence we fix D = E Tq 1 C(i],S) and, with this choice, the conditions 
imposed to get the above bounds, are all implied by: 

E C(r],S) (E C{n,5)\*l 



r \ B5 ' V TqXqS 

for some small enough c'. Hence fixing an order N the induction works for k < AT 
provided we take 8 = £o/4A/" and we find (since AEqC(t}, <5)/(r Ao£o) > 1, for n small 
enough) : 

\ to^o4o ' J- o 

The condition on e has to be such that the map (3.2) can be defined in the domain 

^Co/2,poe- 2 eo,K e- 2 «0' which means that: 

«£ 0(0 o, |c^$|«£o, |^*| <4&, |a,*|«4& (A1.17) 

Bounding the sum $ = £ J$W) and £ j G 0) via the bound (Al.l), (A1.16) 

this means that the corresponding conditions for the existence of the canonical 
maps become: \eBAf\ < 1/2, \sD\ < £, 2 \q, which requires e to be smaller than 
the minimum e between (ToXo^E^C^, S)~ 1 ) 2 Af~ 3 b , with b suitably small, and 

where one power of Af in the factor A/" -3 comes from A/1 < Af M (see (Al.l)) and the 
other two from the Af in (A1.16); and we use also £o < 1 an d r AoC(?7, S)~ 1 E^ 1 < 1 
(see (2.2)), so that the minimum in (A1.18) is reached in the first term. 
It remains to study the remainders of order > TV in the Taylor expansion: this 
can easily done in terms of A", B and D, and a bound 2(e / £o)^ Eq is obtained. 
Then from the fact that C{n, 5) < coTT^o^C^o" 1 ) 2 M-^ 1 ), for M < C^A)" 1 , 
theorem 1 follows, with b' = boc^ 2 small enough. 

A corollary of theorem 1 is a fast averaging result: namely if there are no slow 
frequencies (i.e. if the Fourier transform in a contains only Fourier modes which 
are non-zero multiples of a given u = (0,1/2) or more generally of js = (y\,vi) 
with vi ^ 0, i.e. if / is "monochromatic") the lower bound on A is T C(r], S)" 1 with 
C(t],5) = (772 S~ 2 + 5~ 1 )log5~ 1 . Hence if Af is fixed arbitrarily but rj -independent 
the condition for casting the Hamiltonian in the form (4.4) is |e| < 0(1). 
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